Systematic analysis of RNA-binding proteins identifies targetable therapeutic vulnerabilities in osteosarcoma.

Osteosarcoma is the most common primary malignant bone tumor with a strong tendency to metastasize, limiting the prognosis of affected patients. Genomic, epigenomic and transcriptomic analyses have demonstrated the exquisite molecular complexity of this tumor, but have not sufficiently defined the underlying mechanisms or identified promising therapeutic targets. To systematically explore RNA-protein interactions relevant to OS, we define the RNA interactomes together with the full proteome and the transcriptome of cells from five malignant bone tumors (four osteosarcomata and one malignant giant cell tumor of the bone) and from normal mesenchymal stem cells and osteoblasts. These analyses uncover both systematic changes of the RNA-binding activities of defined RNA-binding proteins common to all osteosarcomata and individual alterations that are observed in only a subset of tumors. Functional analyses reveal a particular vulnerability of these tumors to translation inhibition and a positive feedback loop involving the RBP IGF2BP3 and the transcription factor Myc which affects cellular translation and OS cell viability. Our results thus provide insight into potentially clinically relevant RNA-binding protein-dependent mechanisms of osteosarcoma.

Improved discovery of RNA-binding protein binding sites in eCLIP data using DEWSeq.

Enhanced crosslinking and immunoprecipitation (eCLIP) sequencing is a method for transcriptome-wide detection of binding sites of RNA-binding proteins (RBPs). However, identified crosslink sites can deviate from experimentally established functional elements of even well-studied RBPs. Current peak-calling strategies result in low replication and high false positive rates. Here, we present the R/Bioconductor package DEWSeq that makes use of replicate information and size-matched input controls. We benchmarked DEWSeq on 107 RBPs for which both eCLIP data and RNA sequence motifs are available and were able to more than double the number of motif-containing binding regions relative to standard eCLIP processing. The improvement not only relates to the number of binding sites (3.1-fold with known motifs for RBFOX2), but also their subcellular localization (1.9-fold of mitochondrial genes for FASTKD2) and structural targets (2.2-fold increase of stem-loop regions for SLBP. On several orthogonal CLIP-seq datasets, DEWSeq recovers a larger number of motif-containing binding sites (3.3-fold). DEWSeq is a well-documented R/Bioconductor package, scalable to adequate numbers of replicates, and tends to substantially increase the proportion and total number of RBP binding sites containing biologically relevant features.

Global gene regulatory network underlying miR165a in Arabidopsis shoot apical meristem.

Arabidopsis microRNA165a (miR165a) targets Class III Homeodomain Leucine-Zipper (HD-ZIPIII) transcription factors to regulate various aspects of plant development and stress response. Over-expression of miR165a mimics the loss-of-function phenotype of HD-ZIPIII genes and leading to ectopic organ formation, shoot apical meristem (SAM) termination, loss of leaf polarity, and defective vasculature development. However, the molecular mechanisms underlying these phenotypes remain unresolved. Here, we over-expressed miR165a in a dexamethasone inducible manner and identified differentially expressed genes in the SAM through RNA-Seq. Simultaneously, using multi-channel FACS combined with RNA-Seq approach, we characterized global transcriptome patterns in miR165a expressing cell-types compared to HD-ZIPIII expressing cell-types and other cell-types in SAM. By integrating our results we identified sets of genes which are up-regulated by miR165a as well have enriched expression in miR165a cell-types, and vice-versa. Known plant development related genes such as HD-ZIPIII and their targets LITTLE ZIPPERs, Like AUXIN RESISTANT 2, BEL1-like homeodomain 6, ROTUNDIFOLIA like 16 were found to be down-regulated. Among the up-regulated genes, GIBBERELLIN 2-OXIDASEs, various elemental transporters (YSL3, ZIFL1, SULTR), and other transporter genes were prominent. Thus, the genes identified in this study help to unravel the molecular mechanism of miR165a and HD-ZIPIII regulated plant development and stress response.

The RNA-binding protein landscapes differ between mammalian organs and cultured cells.

System-wide approaches have unveiled an unexpected breadth of the RNA-bound proteomes of cultured cells. Corresponding information regarding RNA-binding proteins (RBPs) of mammalian organs is still missing, largely due to technical challenges. Here, we describe ex vivo enhanced RNA interactome capture (eRIC) to characterize the RNA-bound proteomes of three different mouse organs. The resulting organ atlases encompass more than 1300 RBPs active in brain, kidney or liver. Nearly a quarter (291) of these had formerly not been identified in cultured cells, with more than 100 being metabolic enzymes. Remarkably, RBP activity differs between organs independent of RBP abundance, suggesting organ-specific levels of control. Similarly, we identify systematic differences in RNA binding between animal organs and cultured cells. The pervasive RNA binding of enzymes of intermediary metabolism in organs points to tightly knit connections between gene expression and metabolism, and displays a particular enrichment for enzymes that use nucleotide cofactors. We describe a generically applicable refinement of the eRIC technology and provide an instructive resource of RBPs active in intact mammalian organs, including the brain.

Small noncoding RNA interactome capture reveals pervasive, carbon source-dependent tRNA engagement of yeast glycolytic enzymes.

Small noncoding RNAs fulfill key functions in cellular and organismal biology, typically working in concert with RNA-binding proteins (RBPs). While proteome-wide methodologies have enormously expanded the repertoire of known RBPs, these methods do not distinguish RBPs binding to small noncoding RNAs from the rest. To specifically identify this relevant subclass of RBPs, we developed small noncoding RNA interactome capture (snRIC2C) based on the differential RNA-binding capacity of silica matrices (2C). We define the S. cerevisiae proteome of nearly 300 proteins that specifically binds to RNAs smaller than 200 nt in length (snRBPs), identifying informative distinctions from the total RNA-binding proteome determined in parallel. Strikingly, the snRBPs include most glycolytic enzymes from yeast. With further methodological developments using silica matrices, 12 tRNAs were identified as specific binders of the glycolytic enzyme GAPDH. We show that tRNA engagement of GAPDH is carbon source-dependent and regulated by the RNA polymerase III repressor Maf1, suggesting a regulatory interaction between glycolysis and RNA polymerase III activity. We conclude that snRIC2C and other 2C-derived methods greatly facilitate the study of RBPs, revealing previously unrecognized interactions.

htseq-clip: a toolset for the preprocessing of eCLIP/iCLIP datasets.

SUMMARY: Transcriptome-wide detection of binding sites of RNA-binding proteins is achieved using Individual-nucleotide crosslinking and immunoprecipitation (iCLIP) and its derivative enhanced CLIP (eCLIP) sequencing methods. Here, we introduce htseq-clip, a python package developed for preprocessing, extracting and summarizing crosslink site counts from i/eCLIP experimental data. The package delivers crosslink site count matrices along with other metrics, which can be directly used for filtering and downstream analyses such as the identification of differential binding sites. AVAILABILITY AND IMPLEMENTATION: The Python package htseq-clip is available via pypi (python package index), bioconda and the Galaxy Tool Shed under the open source MIT License. The code is hosted at https://github.com/EMBL-Hentze-group/htseq-clip and documentation is available under https://htseq-clip.readthedocs.io/en/latest.

A Pipeline for Analyzing eCLIP and iCLIP Data with Htseq-clip and DEWSeq.

Individual-nucleotide crosslinking and immunoprecipitation (iCLIP) sequencing and its derivative enhanced CLIP (eCLIP) sequencing are methods for the transcriptome-wide detection of binding sites of RNA-binding proteins (RBPs). This chapter provides a stepwise tutorial for analyzing iCLIP and eCLIP data with replicates and size-matched input (SMI) controls after read alignment using our open-source tools htseq-clip and DEWSeq. This includes the preparation of gene annotation, extraction, and preprocessing of truncation sites and the detection of significantly enriched binding sites using a sliding window based approach suitable for different binding modes of RBPs.

MPP8 is essential for sustaining self-renewal of ground-state pluripotent stem cells.

Deciphering the mechanisms that control the pluripotent ground state is key for understanding embryonic development. Nonetheless, the epigenetic regulation of ground-state mouse embryonic stem cells (mESCs) is not fully understood. Here, we identify the epigenetic protein MPP8 as being essential for ground-state pluripotency. Its depletion leads to cell cycle arrest and spontaneous differentiation. MPP8 has been suggested to repress LINE1 elements by recruiting the human silencing hub (HUSH) complex to H3K9me3-rich regions. Unexpectedly, we find that LINE1 elements are efficiently repressed by MPP8 lacking the chromodomain, while the unannotated C-terminus is essential for its function. Moreover, we show that SETDB1 recruits MPP8 to its genomic target loci, whereas transcriptional repression of LINE1 elements is maintained without retaining H3K9me3 levels. Taken together, our findings demonstrate that MPP8 protects the DNA-hypomethylated pluripotent ground state through its association with the HUSH core complex, however, independently of detectable chromatin binding and maintenance of H3K9me3.

An integrated analysis of cell-type specific gene expression reveals genes regulated by REVOLUTA and KANADI1 in the Arabidopsis shoot apical meristem.

In the Arabidopsis thaliana shoot apical meristem (SAM) the expression domains of Class III Homeodomain Leucine Zipper (HD-ZIPIII) and KANADI (KAN) genes are separated by a narrow boundary region from which new organs are initiated. Disruption of this boundary through either loss of function or ectopic expression of HD-ZIPIII and KAN causes ectopic or suppression of organ formation respectively, raising the question of how these transcription factors regulate organogenesis at a molecular level. In this study we develop a multi-channel FACS/RNA-seq approach to characterize global patterns of gene expression across the HD-ZIPIII-KAN1 SAM boundary. We then combine FACS, RNA-seq and perturbations of HD-ZIPIII and KAN expression to identify genes that are both responsive to REV and KAN1 and normally expressed in patterns that correlate with REV and KAN1. Our data reveal that a significant number of genes responsive to REV are regulated in opposite ways depending on time after induction, with genes associated with auxin response and synthesis upregulated initially, but later repressed. We also characterize the cell type specific expression patterns of auxin responsive genes and identify a set of genes involved in organogenesis repressed by both REV and KAN1.

Transcriptome profile of lung dendritic cells after in vitro porcine reproductive and respiratory syndrome virus (PRRSV) infection.

The porcine reproductive and respiratory syndrome (PRRS) is an infectious disease that leads to high financial and production losses in the global swine industry. The pathogenesis of this disease is dependent on a multitude of factors, and its control remains problematic. The immune system generally defends against infectious diseases, especially dendritic cells (DCs), which play a crucial role in the activation of the immune response after viral infections. However, the understanding of the immune response and the genetic impact on the immune response to PRRS virus (PRRSV) remains incomplete. In light of this, we investigated the regulation of the host immune response to PRRSV in porcine lung DCs using RNA-sequencing (RNA-Seq). Lung DCs from two different pig breeds (Pietrain and Duroc) were collected before (0 hours) and during various periods of infection (3, 6, 9, 12, and 24 hours post infection (hpi)). RNA-Seq analysis revealed a total of 20,396 predicted porcine genes, which included breed-specific differentially expressed immune genes. Pietrain and Duroc infected lung DCs showed opposite gene expression courses during the first time points post infection. Duroc lung DCs reacted more strongly and distinctly than Pietrain lung DCs during these periods (3, 6, 9, 12 hpi). Additionally, cluster analysis revealed time-dependent co-expressed groups of genes that were involved in immune-relevant pathways. Key clusters and pathways were identified, which help to explain the biological and functional background of lung DCs post PRRSV infection and suggest IL-1ß1 as an important candidate gene. RNA-Seq was also used to characterize the viral replication of PRRSV for each breed. PRRSV was able to infect and to replicate differently in lung DCs between the two mentioned breeds. These results could be useful in investigations on immunity traits in pig breeding and enhancing the health of pigs.

Exploring novel mechanistic insights in Alzheimer's disease by assessing reliability of protein interactions.

Protein interaction networks are widely used in computational biology as a graphical means of representing higher-level systemic functions in a computable form. Although, many algorithms exist that seamlessly collect and measure protein interaction information in network models, they often do not provide novel mechanistic insights using quantitative criteria. Measuring information content and knowledge representation in network models about disease mechanisms becomes crucial particularly when exploring new target candidates in a well-defined functional context of a potential disease mechanism. To this end, we have developed a knowledge-based scoring approach that uses literature-derived protein interaction features to quantify protein interaction confidence. Thereby, we introduce the novel concept of knowledge cliffs, regions of the interaction network where a significant gap between high scoring and low scoring interactions is observed, representing a divide between established and emerging knowledge on disease mechanism. To show the application of this approach, we constructed and assessed reliability of a protein-protein interaction model specific to Alzheimer's disease, which led to screening, and prioritization of four novel protein candidates. Evaluation of the identified candidates showed that two of them are already followed in clinical trials for testing potential AD drugs.

Preliminary study of FMO1, FMO5, CYP21, ESR1, PLIN2 and SULT2A1 as candidate gene for compounds related to boar taint.

An association study between polymorphisms of six genes and boar taint related compounds androstenone, skatole and indole was performed in a boar population (n=370). Significant association (P<0.05) was detected for SNP of FMO5 (g.494A>G) with all boar taint compounds, SNP of CYP21 (g.3911T>C) with skatole and indole, and SNP of ESR1 (g.672C>T) with androstenone and indole. mRNA expression of CYP21 and ESR1 was higher in CAB (castrated boar) compared to non-castrated boars; whereas, the expression of FMO5 and ESR1 was higher in LBT (low boar taint) compared to HBT (high boar taint) in liver tissue. FMO5, CYP21 and ESR1 proteins were less detectable in HBT compared with LBT and CAB in liver tissues. These findings suggest that FMO5, CYP21 and ESR1 gene variants might have effects on the boar taint compounds.

MicroRNA Expression Profile in Bovine Granulosa Cells of Preovulatory Dominant and Subordinate Follicles during the Late Follicular Phase of the Estrous Cycle.

In bovine, ovarian follicles grow in a wave-like fashion with commonly 2 or 3 follicular waves emerging per estrous cycle. The dominant follicle of the follicular wave which coincides with the LH-surge becomes ovulatory, leaving the subordinate follicles to undergo atresia. These physiological processes are controlled by timely and spatially expressed genes and gene products, which in turn are regulated by post-transcriptional regulators. MicroRNAs, a class of short non-coding RNA molecules, are one of the important posttranscriptional regulators of genes associated with various cellular processes. Here we investigated the expression pattern of miRNAs in granulosa cells of bovine preovulatory dominant and subordinate follicles during the late follicular phase of bovine estrous cycle using Illumina miRNA deep sequencing. In addition to 11 putative novel miRNAs, a total of 315 and 323 known miRNAs were detected in preovulatory dominant and subordinate follicles, respectively. Moreover, in comparison with the subordinate follicles, a total of 64 miRNAs were found to be differentially expressed in preovulatory dominant follicles, of which 34 miRNAs including the miR-132 and miR-183 clusters were significantly enriched, and 30 miRNAs including the miR-17-92 cluster, bta-miR-409a and bta-miR-378 were significantly down regulated in preovulatory dominant follicles. In-silico pathway analysis revealed that canonical pathways related to oncogenesis, cell adhesion, cell proliferation, apoptosis and metabolism were significantly enriched by the predicted target genes of differentially expressed miRNAs. Furthermore, Luciferase reporter assay analysis showed that one of the differentially regulated miRNAs, the miR-183 cluster miRNAs, were validated to target the 3'-UTR of FOXO1 gene. Moreover FOXO1 was highly enriched in granulosa cells of subordinate follicles in comparison with the preovulatory dominant follicles demonstrating reciprocal expression pattern with miR-183 cluster miRNAs. In conclusion, the presence of distinct sets of miRNAs in granulosa cells of preovulatory dominant and subordinate follicles supports the potential role of miRNAs in post-transcriptional regulation of genes involved in bovine follicular development during the late follicular phase of the estrous cycle.

Identification of gene co-expression clusters in liver tissues from multiple porcine populations with high and low backfat androstenone phenotype.

BACKGROUND: Boar taint is principally caused by accumulation of androstenone and skatole in adipose tissues. Studies have shown high heritability estimates for androstenone whereas skatole production is mainly dependent on nutritional factors. Androstenone is a lipophilic steroid mainly metabolized in liver. Majority of the studies on hepatic androstenone metabolism focus only on a single breed and very few studies account for population similarities/differences in gene expression patterns. In this work, we concentrated on population similarities in gene expression to identify the common genes involved in hepatic androstenone metabolism of multiple pig populations. Based on androstenone measurements, publicly available gene expression datasets from three porcine populations were compiled into either low or high androstenone dataset. Gene expression correlation coefficients from these datasets were converted to rank ratios and joint probabilities of these rank ratios were used to generate dataset specific co-expression clusters. Finally, these networks were clustered using a graph clustering technique. RESULTS: Cluster analysis identified a number of statistically significant co-expression clusters in the dataset. Further enrichment analysis of these clusters showed that one of the clusters from low androstenone dataset was highly enriched for xenobiotic, drug, cholesterol and lipid metabolism and cytochrome P450 associated metabolism of drugs and xenobiotics. Literature references revealed that a number of genes in this cluster were involved in phase I and phase II metabolism. Physical and functional similarity assessment showed that the members of this cluster were dispersed across multiple clusters in high androstenone dataset, possibly indicating a weak co-expression of these genes in high androstenone dataset. CONCLUSIONS: Based on these results we hypothesize that majority of the genes in this cluster forms a signature co-expression cluster in low androstenone dataset in our experiment and that majority of the members of this cluster might be responsible for hepatic androstenone metabolism across all the three populations used in our study. We propose these results as a background work towards understanding breed similarities in hepatic androstenone metabolism. Additional large scale experiments using data from multiple porcine breeds are necessary to validate these findings.

The expression pattern of microRNAs in granulosa cells of subordinate and dominant follicles during the early luteal phase of the bovine estrous cycle.

This study aimed to investigate the miRNA expression patterns in granulosa cells of subordinate (SF) and dominant follicle (DF) during the early luteal phase of the bovine estrous cycle. For this, miRNA enriched total RNA isolated from granulosa cells of SF and DF obtained from heifers slaughtered at day 3 and day 7 of the estrous cycle was used for miRNAs deep sequencing. The results revealed that including 17 candidate novel miRNAs, several known miRNAs (n = 291-318) were detected in SF and DF at days 3 and 7 of the estrous cycle of which 244 miRNAs were common to all follicle groups. The let-7 families, bta-miR-10b, bta-miR-26a, bta-miR-99b and bta-miR-27b were among abundantly expressed miRNAs in both SF and DF at both days of the estrous cycle. Further analysis revealed that the expression patterns of 16 miRNAs including bta-miR-449a, bta-miR-449c and bta-miR-222 were differentially expressed between the granulosa cells of SF and DF at day 3 of the estrous cycle. However, at day 7 of the estrous cycle, 108 miRNAs including bta-miR-409a, bta-miR-383 and bta-miR-184 were differentially expressed between the two groups of granulosa cell revealing the presence of distinct miRNA expression profile changes between the two follicular stages at day 7 than day 3 of the estrous cycle. In addition, unlike the SF, marked temporal miRNA expression dynamics was observed in DF groups between day 3 and 7 of the estrous cycle. Target gene prediction and pathway analysis revealed that major signaling associated with follicular development including Wnt signaling, TGF-beta signaling, oocyte meiosis and GnRH signaling were affected by differentially expressed miRNAs. Thus, this study highlights the miRNA expression patterns of granulosa cells in subordinate and dominant follicles that could be associated with follicular recruitment, selection and dominance during the early luteal phase of the bovine estrous cycle.

Pathway based analysis of genes and interactions influencing porcine testis samples from boars with divergent androstenone content in back fat.

One of the primary factors contributing to boar taint is the level of androstenone in porcine adipose tissues. A majority of the studies performed to identify candidate biomarkers for the synthesis of androstenone in testis tissues follow a reductionist approach, identifying and studying the effect of biomarkers individually. Although these studies provide detailed information on individual biomarkers, a global picture of changes in metabolic pathways that lead to the difference in androstenone synthesis is still missing. The aim of this work was to identify major pathways and interactions influencing steroid hormone synthesis and androstenone biosynthesis using an integrative approach to provide a bird's eye view of the factors causing difference in steroidogenesis and androstenone biosynthesis. For this purpose, we followed an analysis procedure merging together gene expression data from boars with divergent levels of androstenone and pathway mapping and interaction network retrieved from KEGG database. The interaction networks were weighted with Pearson correlation coefficients calculated from gene expression data and significant interactions and enriched pathways were identified based on these networks. The results show that 1,023 interactions were significant for high and low androstenone animals and that a total of 92 pathways were enriched for significant interactions. Although published articles show that a number of these enriched pathways were activated as a result of downstream signaling of steroid hormones, we speculate that the significant interactions in pathways such as glutathione metabolism, sphingolipid metabolism, fatty acid metabolism and significant interactions in cAMP-PKA/PKC signaling might be the key factors determining the difference in steroidogenesis and androstenone biosynthesis between boars with divergent androstenone levels in our study. The results and assumptions presented in this study are from an in-silico analysis done at the gene expression level and further laboratory experiments at genomic, proteomic or metabolomic level are necessary to validate these findings.

Identification of the novel candidate genes and variants in boar liver tissues with divergent skatole levels using RNA deep sequencing.

Boar taint is the unpleasant odour of meat derived from non-castrated male pigs, caused by the accumulation of androstenone and skatole in fat. Skatole is a tryptophan metabolite produced by intestinal bacteria in gut and catabolised in liver. Since boar taint affects consumer's preference, the aim of this study was to perform transcriptome profiling in liver of boars with divergent skatole levels in backfat by using RNA-Seq. The total number of reads produced for each liver sample ranged from 11.8 to 39.0 million. Approximately 448 genes were differentially regulated (p-adjusted <0.05). Among them, 383 genes were up-regulated in higher skatole group and 65 were down-regulated (p<0.01, FC>1.5). Differentially regulated genes in the high skatole liver samples were enriched in metabolic processes such as small molecule biochemistry, protein synthesis, lipid and amino acid metabolism. Pathway analysis identified the remodeling of epithelial adherens junction and TCA cycle as the most dominant pathways which may play important roles in skatole metabolism. Differential gene expression analysis identified candidate genes in ATP synthesis, cytochrome P450, keratin, phosphoglucomutase, isocitrate dehydrogenase and solute carrier family. Additionally, polymorphism and association analysis revealed that mutations in ATP5B, KRT8, PGM1, SLC22A7 and IDH1 genes could be potential markers for skatole levels in boars. Furthermore, expression analysis of exon usage of three genes (ATP5B, KRT8 and PGM1) revealed significant differential expression of exons of these genes in different skatole levels. These polymorphisms and exon expression differences may have impacts on the gene activity ultimately leading to skatole variation and could be used as genetic marker for boar taint related traits. However, further validation is required to confirm the effect of these genetic markers in other pig populations in order to be used in genomic selection against boar taint in pig breeding programs.

RNA deep sequencing reveals novel candidate genes and polymorphisms in boar testis and liver tissues with divergent androstenone levels.

Boar taint is an unpleasant smell and taste of pork meat derived from some entire male pigs. The main causes of boar taint are the two compounds androstenone (5α-androst-16-en-3-one) and skatole (3-methylindole). It is crucial to understand the genetic mechanism of boar taint to select pigs for lower androstenone levels and thus reduce boar taint. The aim of the present study was to investigate transcriptome differences in boar testis and liver tissues with divergent androstenone levels using RNA deep sequencing (RNA-Seq). The total number of reads produced for each testis and liver sample ranged from 13,221,550 to 33,206,723 and 12,755,487 to 46,050,468, respectively. In testis samples 46 genes were differentially regulated whereas 25 genes showed differential expression in the liver. The fold change values ranged from -4.68 to 2.90 in testis samples and -2.86 to 3.89 in liver samples. Differentially regulated genes in high androstenone testis and liver samples were enriched in metabolic processes such as lipid metabolism, small molecule biochemistry and molecular transport. This study provides evidence for transcriptome profile and gene polymorphisms of boars with divergent androstenone level using RNA-Seq technology. Digital gene expression analysis identified candidate genes in flavin monooxygenease family, cytochrome P450 family and hydroxysteroid dehydrogenase family. Moreover, polymorphism and association analysis revealed mutation in IRG6, MX1, IFIT2, CYP7A1, FMO5 and KRT18 genes could be potential candidate markers for androstenone levels in boars. Further studies are required for proving the role of candidate genes to be used in genomic selection against boar taint in pig breeding programs.